NAOJ-Th-Ap2003, No.22 

Preprint typeset using L^TgK style emulateapj v. 1 1/12/01 



DIRECT MAPPING OF MASSIVE COMPACT OBJECTS IN EXTRAGALACTIC DARK HALOS 
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ABSTRACT 

A significant fraction of non-baryonic or baryonic dark matter in galactic halos may consist of MASsive Com- 
pact Objects (MASCOs) with mass M = IQ^'^Mq. Possible candidates for such compact objects include pri- 
mordial black holes or remnants of primordial (Population III) stars. We propose a method for directly detect- 
ing MASCOs in extragalactic halos, using the VLBI techniques with extremely high resolution that would be 
achieved by the next generation mission of the VLBI Space Observatory Program such as VSOP-2. If a galactic 
halo comprising a large number of MASCOs produces multiple images of a background radio-loud QSO by grav- 
itational lensing, then a high-resolution radio map of each macro-lensed image should reveal microlensing effects 
by MASCOs. To assess their observational feasibility, we simulate microlensing of the radio-loud, four-image 
lensed QSO, B1422-I-231, assuming angular resolution of ^.01 mas. MASCOs are represented by point masses. 
For comparison, we also simulate microlensing of B 1422-1-231 by singular isothermal spheres. We find that the 
surface brightness of the macro-lensed images shows distinct spatial patterns on the scale of the Einstein radius 
of the perturbers. In the case of point-mass perturbers, many tiny dark spots also appear in the macro-lensed 
images associated with a decrease in the surface brightness toward the fringe of the original QSO image, whereas 
no such spots are available in the SIS models. Because such spatial patterns in each macro-lensed image cannot 
be linearly mapped to those in other macro-lensed images if they are relevant to lensing perturbers, it is fairly 
easy to discriminate them from intrinsic substructures within a QSO. Based on the size, position and magnified 
or demagnified patterns of images, we shall be able to determine the mass and density profile of an individual 
MASCO as well as its spatial distribution and abundance in a galactic halo. 

Subject headings: gravitational lensing - dark matter - galaxies: cluster - stars: Population III 



1. INTRODUCTION 

Gravitational lensing is a powerful tool for probing dark mat- 
ter in the form of compact objects, such as black holes, neu- 
tron stars, and dwarfs. To date, various kinds of gravitational 
lensing studies give stringent upper limits on the abundance of 
these objects on very different mass scales, ranging from sub- 
stellar objects 1O~^M0 (Alcock et al. 2000) to galaxy clusters 
- IO'^^Mq (Nemiroff 1991). 

However, there have been no stringent observational con- 
straints on dark compact objects with a mass range of 1O'M0 < 
M < 1O'*M0 corresponding to the mass range of massive stars. 
For brevity, we call them MASCOs (MASsive Compact Ob- 
jects). It is of particular importance to constrain their abun- 
dance in galactic dark halos, because MASCOs are both bary- 
onic and non-baryonic dark matter candidates. Recent obser- 
vation of the cosmic microwave background (Bennett et al. 
2003) suggests that the primordial stars (or Population III stars) 
are rather massive as indicated by theoretical works (Bromm, 
Coppi, & Larson 1999; Abel, Bryan, & Norman 2000; Tsuribe 
& Inutsuka 2001 ; Omukai & Palla 2002) and their formation ef- 
ficiency are higher than current theoretical works suggest (Cen 
2003). In some MASCO mass scales, primordial stars collapse 
directly to form black holes (Bond et al. 1984; Woosley 1986). 
Thus, a large fraction of baryonic dark matter in galactic dark 
halos may consist of black holes that are remnants of massive 
primordial stars. Alternatively, as the cold dark matter (CDM) 
candidates, MASCOs may be primordial black holes that are 
produced in the early universe at the era between the QCD 
phase transition and the big-bang nucleosynthesis (Carr 1975). 

To date, various gravitational lensing methods that can probe 
MASCOs have been proposed. Gould (1992) proposed that an- 



nual oscillation in the light curve of a MACHO lensing event 
can distinguish a MASCO from a stellar mass compact object. 
However, for a mass range of lO^M© <M< IO^Mq, the annual 
oscillation owing to the Earth's motion is too small to determine 
the transverse velocity. Totani (2003) argued that observation 
of cluster-cluster microlensing system by a 8 m-class ground- 
based telescope can probe the abundance of compact objects 
with mass lO-^M© <M < IO'^Mq. Although his proposed 
method covers a wide range of mass scales, the method itself 
cannot break the degeneracy in the lensing parameters such as 
mass or transverse velocity owing to the lack of information of 
the source star luminosity. 

As an alternative, we propose a method to directly detect 
MASCOs using the VLBI techniques with extremely high res- 
olution. Suppose that a distant radio-loud QSO shows multiple 
images owing to the lensing effect of a foreground galaxy. If 
there are a large number of MASCOs in the galactic halo, then 
the high-resolution radio map of the gravitationally lensed QSO 
images should reveal the microlensing effects by MASCOs. In 
contrast to the previous methods, our proposed method can de- 
termine the mass of MASCOs in the dark halo without any as- 
sumption on the lensing parameters. Furthermore, we can also 
prove the spatial distribution of MASCOs in the dark halo pro- 
jected onto the surface perpendicular to the line of sight. A pos- 
sible target lens system is B 1422+231 in which the source is a 
radio-loud QSO at redshift Zs = 3.62 and the lens is an elliptical 
galaxy at redshift zl = 0.34. Although the method we propose 
here is similar to that of Wambsganss & Paczy'nski (1992) in 
which the target compact object mass is set to M ^ 10''Mq, the 
lensing effect by MASCOs can be quite different since the Ein- 
stein radius is much smaller than the size of a source and the 
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total number in a dark halo is large. In this Letter, we show our 
simulation results for the lens system B 1422+231 microlensed 
by MASCOs, that would be observable by the future VLBI net- 
work. 

Wc propose to construct a new VLBI network with ex- 
tremely high resolution such as the planned next generation 
VLBI Space Observatory Program (VSOP-2) (Hirabayashi et 
al. 2001) and produce high resolution maps of strongly lensed 
radio-loud QSOs such as B 1422+231. Such maps might pro- 
vide us with valuable information of the first star remnants and 
the initial condition of density fluctuations on MASCO scales 
as well as the nature of the CDM. 

2. METHOD 

We consider a lensing system in which a radio-loud QSO is 
strongly lensed by the dark halo of an intervening galaxy. If 
MASCOs constitute a large amount of the dark matter compo- 
nent, then we expect to see many microlensing events by them, 
because the microlensing optical depth is nearly 1 in a strong 
lens. In order to detect the signatures of such microlensing 
in the macro-lensed QSO images, one needs extremely high- 
resolution radio mapping of the images. The necessary resolu- 
tion can be estimated by the Einstein radius of MASCOs 

where Di, Ds, and Dis are the angular diameter distance to 
the lens, source, and distance to the lens from the source, re- 
spectively. Although current observations cannot achieve such 
a high-resolution, it is not out of reach in the near future. 
For instance, the target resolution of the VSOP-2 mission is 
^ 0.025 mas at 43 GHz, which can resolve MASCOs with mass 
M > IO^Mq for lensing systems at 1 Gpc. The VSOP-2 satel- 
lite will be placed at an elUptical orbit with an apogee height of 
~30,000 km and a perigee height of ~ 1,000 km (Hirabayashi 
et al. 2001). The resolution can be further improved by a factor 
of two or three by placing satellite(s) at a much distant orbit or 
by increasing the observational frequency. 

In order to demonstrate the microlensing effects of MAS- 
COs, we consider the QSO lensing system B 1422+231 whose 
images consist of three highly magnified ones A, B, and C, and 
a faint one D located near the lens galaxy. The redshifts of the 
source and the lens are zs = 3.62 and zl = 0.34, respectively. In 
what follows, we assume the following cosmological parame- 
ters: ri„, = 0.3, Oa = 0.7 and h = 0.7, giving = 1.00 Gpc and 
£)s= 1-49 Gpc. 

To model the macro-lensing system, we adopt a singular 
isothermal elUpsoid (SIE) in an external shear field in which 
the isopotential curves in the projected surface perpendicular 
to the line of sight are ellipse (Kormann et al. 1994). This 
lens model has the following parameters: Einstein angular ra- 
dius 6e, a ratio q of the minor axis to the major axis of the 
eUipse of the iso-potential contour, position angle (pg of the mi- 
nor axis, source position on the source plane, lens position on 
the lens plane, unlensed source flux, strength and direction of 
the external shear ij,4>s)- The best-fit lensing parameters us- 
ing fit is Oe = 0.78 arcsec, q = 0.84, (j)g = -57.5°, 7 = 0.21, 
(ps = -53.9° (Chiba 2002), where the one-dimensional velocity 
dispersion of the lens is estimated as cr = 183.23 km s~'. The 
observed position of each image is taken from HST observa- 
tions with the Faint Object Camera by Impey et al (1996). We 
show, in the middle panel of Figure 1, the model results for the 



four lensed images. It should be noted that B 1422+231 shows 
an anomalous flux ratio in the macro-lensed images possibly 
caused by the CDM subhalos with masses larger than IO^Mq 
(Chiba 2002). We neglect below the effects of such subhalos 
which are nol relevant to microlensing by MASCOs. 

We estimate the surface number density of MASCOs with 
an individual mass of M in the lens galaxy. Because the ellip- 
ticity of the lens in B 1422+231 is modest, it can be approxi- 
mated by that of a singular isothermal sphere (SIS) with sur- 
face mass density S(/?) = a^/2GR where R is the projected 
radial distance from the lens center. The radial distances of 
the images A, B, and C from the center are roughly equal to 
the linear Einstein radius, ^0 = D^Oe = 3.8 kpc. Note that 
1 mas corresponds to ^4.9 pc on the lens plane. Then the 
surface number density of MASCOs at = ^0 is estimated as 
NiCo) = /S(^o) /M ~ / X 2.4 X m\M/ IO^Mq)-^ mas'^, where 
/ < 1 is a parameter that represents the fraction of MASCOs 
in the total mass of the dark halo. For the D image, the surface 
number density will be 2.7 times larger because the radial dis- 
tance is as small as ^ = £,o/2.1. Assuming that the observable 
radio-emitting region of the QSO has a typical Unear size of the 
order of ~ 10 pc, then the total number of MASCOs overlaid 
on the macro-lensed images would be several for M = IO^Mq 
even if the MASCO fraction is of the order of a few percents. 
The presence of more extended jets in the images will increase 
this number. 

First, we assume that MASCOs are placed randomly on the 
lens plane in which the mean number is determined by the sur- 
face number density N{R). Clustering of MASCOs on sub- 
galactic scales will not be considered for the sake of simpUcity; 
if a clustering signature is observed in the lensed images, it will 
give us a valuable information about the formation process of 
MASCOs. 

Second, we calculate the total gravitational potential of the 
lens from the sum of a SIE that best fits the position and the flux 
ratios of the four macro-lensed images and perturbers, whereas 
the consistent calculation should arrange this combined lens po- 
tential rather than the potential of the SIE alone. However, 
this simpUfication does not affect the concerned microlensing 
signatures by MASCOs. In our simulation, MASCOs are rep- 
resented by point masses which perturb the macro-lens. We 
also carry out simular calculation in which the perturbers are 
represented by SISs for comparison. We also assume that the 
unlensed QSO image has a circular symmetry with a radial pro- 
file given by the Gaussian distribution with a standard deviation 
3ri pc. 

3. RESULTS 

Four panels in Figure 1 show the portions of four macro im- 
ages disturbed by the point-mass MASCOs with M = 10^ Mq 
and / = 0.06, where a resolution of 0.01 mas is assumed. One 
can clearly see several dark spots in the brightest part of the 
macro images as well as several neighboring pairs of dark and 
bright spots in the region where the surface brightness of the 
image has a spatial gradient. These characteristic properties of 
microlensing by point masses can be explained in the following 
manner. 

Let a point mass be located at x = on the lens plane for the 
extended image on the source plane with coordinate y, where 
both X and y are measured in units of the Einstein radius of a 
MASCO. The lens equation is then given by y = x-x/x^. It 
follows that the positions on the lens plane within the Einstein 
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Fig. 1. — A simulation of radio maps of B1422+231 smoothed by a Gaussian beam with FWHIVtO. Olmas. MASCOs are represented by point masses with mass 
M = IO^Mq. Dark spots in zoomed four images correspond to MASCOs in the dark halo of the foreground galaxy. The iso-surface brightness contour interval is 
1/40 of the maximum value. The MASCO fraction to the dark matter halo is assumed to be /=0.06 corresponding to MASCO density ~0.018 provided that the 
total matter density is =0.3. 




Fig. 2. — Simulations of radio maps of B image in B1422+231 smoothed by a Gaussian beam with FWHM=0. Olmas in which singular isothermal spheres with 
mass M = IO'^Mq are placed randomly in the dark halo. The ratio of the SIS radius to the tidal disruption radius is l(left) and l/4(right). 
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radius satisfying -0.62 < x < are mapped to y > 1 . Con- 
versely, the source image extended outside the Einstein radius 
is mapped to the lens-plane region with radius x = 0.62 inside 
the Einstein radius. Suppose that the surface brightness of the 
source has a profile (7,, = y''^ with /? > 0. Then, in the neigh- 
borhood of the lens center x = e <C 1 with e > 0, the surface 
brightness of the lensed image, a^, is decreased as ax = e^. 
Thus, if Uy decreases monotonically with increasing distance 
from the lens center, then one would observe a dark spot (inside 
the Einstein radius) at the lens center. If the sign of is minus, 
then one would observe a bright spot within the Einstein radius. 
Therefore, there will be a neighboring pair of dark and bright 
spots (within the Einstein radius) in the region mapped from 
the source plane where the gradient of ay has a non-vanishing 
constant value. For less massive MASCOs, such discontinuous 
changes of the lensed image become inconspicuous because the 
change in the surface brightness on the scale of the Einstein ra- 
dius becomes smaller. However, if the source has an intrinsic 
substructure, it may have a locally large gradient in its surface 
brightness, yielding pairs of dark and bright spots in the lensed 
image even for less massive MASCOs. It is worth noting that 
we can easily distinguish these microlensing effects from any 
intrinsic structures of the source by linearly mapping the macro- 
lensed image to one of the other images. 

Figure 2 shows the microlensing effects of the SISs with 
M = W^Mq on the B image. In this case, the size of each SIS, 
rsis, is a free parameter, affecting the strength of the lensing ef- 
fect. Its upper limit is set by a tidal radius, r, = (GMr^ /na^y^^, 
at its position r from the galaxy center. We examine two cases, 
''SIS /'"f = 1 and 1 /4, when r is defined as the projected distance 
of the image B from the galaxy center, i.e., r = ^o- In the SIS 
model, the lens equation is given as y = x-x/x. In sharp con- 
trast to the point mass model, the source image extended out- 
side the Einstein radius y > I cannot be mapped to the region 
inside the Einstein radius x < 1, indicating that there will be 
almost no visible effect from a lens placed in the region where 
the gradient in the surface brightness within the Einstein radius 
vanishes. On the other hand, in the region where the surface 
brightness gradient does not vanish, a discontinuous change of 
the image brightness on the scale of the Einstein radius is ob- 
served, as found for the point-mass model. However, magni- 
fication or demagnification effect by SISs is less conspicuous 
than the point-mass model. 

4. SUMMARY AND DISCUSSION 

In this Letter, we proposed a method to directly ditect MAS- 
COs, i.e., massive compact objects with mass M = lO'^M© 



in extragalactic dark halos. Using the VLBl techniques, radio 
maps of strongly lensed radio-loud QSO with extremely high 
resolution ~ 0.01 mas would reveal the abundance, the spatial 
distribution, and the mass density profile of the MASCOs. 

We present numerical simulations for B1422H-231, a typical 
strongly lensed QSO system. In the presence of MASCOs, tiny 
dark spots and pairs of dark and bright spots will appear in the 
macro-lensed images. From the size of the spot, one can esti- 
mate the mass of the MASCO. Even if the total mass of MAS- 
COs in the dark halo of the lens galaxy constitutes only a few 
percents of the total dark halo mass, we can detect MASCOs 
from characteristic substructures in the map of the lensed QSO 
image which cannot be linearly mapped to other macro-lensed 
images. If a perturbing object of the macro-lensed image is de- 
scribed by a SIS, the signature of microlensing appears only if 
the gradient in the surface brightness is very large or the size 
of the SIS is sufficiently compact. In contrast to the point-mass 
case, light from the fringe of the source is not visible within the 
Einstein radius. Therefore, one can distinguish SISs from point 
masses by analyzing the detailed spatial variation in the surface 
brightness. 

If the WMAP result of the Thompson optical depth with- 
stands future data, we will have strong observational evidence 
that the primordial stars are massive and their formation rate 
is quite large (Cen 2003). For mass ranges of SOM© < M < 
IOOM0 and 26OM0 < M, the progenitor of the primordial gas 
can directly collapse to form massive black holes without eject- 
ing any metals (Woosley 1986). Thus the primordial star rem- 
nants may constitute a large fraction of the baryonic dark mat- 
ter. Because our proposed method is so powerful in constrain- 
ing the abundance of point-mass MASCOs, we might be able 
to constrain or determine the initial mass function of the pri- 
mordial stars from the observed mass and the abundance of the 
MASCOs. 

If the CDM is strongly clustered on MASCO mass scales, 
their subhalos might be detected by our proposed method if 
they are compact enough. In our simulation, with resolution 
of 0.01 mas, SIS with mass M ~ 10*Mq can be detected if the 
ratio of the tidal disruption radius to the SIS cutoff radius is 
< 4. By numerically estimating the abundance of subhalos on 
MASCO scales, one can constrain the amplitude of the initial 
density fluctuations on such small scales from future radio ob- 
servations. 

We thank S. Kameno for useful comments on the specifica- 
tion of VSOP-2. This work was supported in part by Grant-in- 
Aid for Scientific Research Fund (No. 1 1367). 
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